function dx=Plant(t,x,flag,para)
dx=zeros(2,1);

ut=para(1);
dt=para(2);


dx(1)=x(2);
dx(2)=-2*x(1)-3*(1-x(1)^2)*x(2)+ut+dt;
